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Abstract 



A recent global analysis of direct photon production at hadron collider 
and fixed target experiments has noted a disturbing trend of disagreement 
between next-to-leading-order (NLO) calculations and data. The conjecture 
has been made that the discrepancy is due to explicit multiple parton emission 
effects which are not accounted for in the theoretical calculations. We investi- 
gate this problem by merging a NLO calculation of direct photon production 
with extra multiple parton emissions via the parton shower (PS) algorithm. 
Our calculation maintains the integrity of the underlying NLO calculation 
while avoiding ambiguities due to double counting of multiple parton emis- 
sions. We find that the NLO+PS calculation can account for much of the 
theory/CDF data discrepancy at y/s = 1.8 TeV. It can also account for much 
of the theory/UA2 discrepancy if a very large virtuality is assumed to initiate 
the initial state parton shower. For lower energy data sets (e.g. y/s < 63 
GeV), NLO+PS calculations alone cannot account for the data/theory dis- 
crepancy, so that some additional non-perturbative kx smearing is needed. 

PACS numbers: 12.38.Bx, 12.38. Qk, 13.85.Qk, 14.80.Am 
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I. INTRODUCTION 



Direct photon production [|I]-|JJ in hadronic collisions has long been recognized as an 
important testing ground for perturbative QCD since many of the ambiguities involved in 
measuring jets are not present when analyzing photons. Direct photon production in lowest 
order QCD takes place via annihilation and Compton scattering Feynman graphs. Since 
the Compton graph involves initial state gluon-quark scattering, measurements of direct 
photon events can serve as important constraints in the determination of the gluon parton 
distribution function ||. For such a program to proceed, the greater precision involved 
in next-to-leading order (NLO) QCD calculations for the hard scattering are used. NLO 
calculations for parton+parton^ 7X have been performed both analytically || and in a 
Monte Carlo framework jTJ. 

A recent global analysis of direct photon production in hadron collisions has noted a 
discrepancy between NLO calculations and a large array of data for the transverse momen- 
tum px distributions of the photon || . Characteristically, in both fixed target and collider 
experiments, there is an experimental excess of photons at low transverse momentum. Sev- 
eral possible explanations have been put forth to resolve the discrepancy. These include i) 
improved (NLO) treatment of bremsstrahlung contributions @ and isolation criteria ||10|| ) ii) 
modifying gluon distribution functions and QCD scale choices to improve the data/theory 
agreement [10|, or usage of alternative parton distribution functions (PDF's) |pT| , and Hi) 



invoking additional partonic kj- smearing effects ||. The latter case comes in two different 
guises: extra partonic kx can come from non-perturbative effects from parton binding and 
intrinsic transverse momentum, or from additional hard multiple parton emissions which 
can be calculated or modeled in perturbative QCD. The non-perturbative effects are gener- 
ally implemented as Gaussian smearing in an attempt to match the data. The perturbative 
multiple gluon emission effects can be implemented via even higher (but fixed) order per- 
turbative calculations, via multiple gluon resummation techniques, or via the parton shower 
(PS) algorithm |fl2|JT3| . The resummation and PS approaches both involve approximate all 
orders perturbative QCD effects. 

In this paper, we explore the extent to which the direct photon data can be explained 
by combining a NLO QCD calculation with multiple parton emission via the parton shower 
algorithm. In doing so, we follow generally the prescription outlined in Ref. |T3j], where NLO 



W and Z boson production were merged with parton showers. In these calculations, Owens' 
method of phase space slicing is used to evaluate the NLO cross sections [|HJ. This method 
lends itself to a direct implementation of parton showers wherein a potential problem of 
double counting multiple parton emissions can be avoided. We show that our implemen- 
tation of showering with the NLO QCD calculation yields an excess of events at low px 
relative to the unshowered NLO result at the Fermilab Tevatron and CERN SppS energies, 
qualitatively accounting for the discrepancy between theory and experiment. Additional 
non-perturbative smearing is required for lower energies characteristic of the CERN ISR or 
fixed target experiments. 
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II. CALCULATIONAL PROCEDURE 



Central to our calculation of direct photon production is the numerical integration of 
phase space via Monte Carlo methods 0. One begins by evaluating the 0(aa s ) and 0(aaf) 
direct photon production subprocess Feynman graphs, including bremsstrahlung corrections 
to qq — > qq, etc. Dimensional regularization is used here for ultraviolet, soft and collinear 
singularities. The four-momenta for the 2 — > 2 subprocess are labeled according to, for 
instance, g{j>\) + q{j>2) — > l(j>z) + i(Pa)] similarly, for 2 — > 3 subprocesses, we use g{p\) + 
l(P2) — > l{Ps) + <?G°4) + 9(1*5)1 e tc- Ultraviolet singularities are renormalized using the 
MS prescription ||. Collinear singularities are factorized and then absorbed into parton 
distribution functions (PDF's) or fragmentation functions. Soft singularities are canceled 
between 2^3 graphs and 2^2 graphs. At this point, all cross section contributions are 
finite, so that numerical predictions can be made. 

What is peculiar to the Monte Carlo method of NLO calculation used here is that the 
phase space integrations are done partly analytically, and partly numerically. The boundary 
between numerical and analytical methods is chosen by selecting two theoretical cutoffs to 
demarcate the collinear and soft regimes. If any invariant quantity tij = (pi — Pj) 2 from 
the 2 — > 3 subprocess has a value < <5 c si2, where Sjj = (pi +pj) 2 , then one is in the 
collinear regime. In this regime, the matrix element squared is evaluated in the leading pole 
approximation and the integration near the collinear pole is done analytically. The cross 
section contribution is de facto 2^2, and it is combined with the leading order and virtual 
contributions to the 2^2 subprocesses. If a final state gluon energy (in the subprocess 
rest frame) has value E g < 5 sy /si2/2, then one is in the soft regime. The integrations of 
the squared matrix elements are performed analytically using the soft gluon approximation, 
and combined with contributions from 2 — > 2 subprocesses. The total 2 — > 2 results, after 
factorization, are finite, but depend on 5 S and 5 C , such that the soft and collinear singularities 
are recovered in the 5 — > limit. The remaining phase space integrations are performed via 
Monte Carlo. This allows easy binning of any desired observables and allows for the simple 
evaluation of the effect of experimental cuts on the NLO prediction ]7j[I3). The 2-^3 
contributions are all positive definite over phase space, but are also singular as 5 S — > or 
5 C — > 0. The 2^2 contributions compensate the 2^3 contributions and result in a 
total cross section which is independent of 5 S and 5 C over a wide range of values M. The 
expressions for all 2 — * 2 and 2 — > 3 processes in direct photon production, through NLO, are 
compiled in Ref . . This is the starting point of our evaluation of the transverse momentum 
of the direct photon using a merger of NLO QCD and parton showers. 

The PS algorithm combines the simplified collinear dynamics, represented by the Q 2 
evolution of parton distribution functions and fragmentation functions, with the exact kine- 



matics of multiple parton emission |l^JT3|] . As implemented here, no additional weights 
to the integral are included with parton showers, as the Q 2 evolved distribution functions 
and fragmentation functions are used in evaluating the differential cross section. For the 
direct photon transverse momentum distribution, initial rather than final state showering 



is most important. Using a backward shower algorithm [fT5 ], the initial state showers are 
evolved backward from a starting virtuality t v . The kinematics of the multiple partons in 
the initial state shower result in transverse momenta for the partons participating in the 
hard scattering, effectively boosting the direct photon transverse momentum relative to the 
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collinear approximation of the kinematics. In practice, the parton shower is cutoff at some 
low t min value where perturbative QCD is still valid, but where the multiple emissions no 
longer become resolvable. In all the results described below, we set t min = 5 GeV 2 . Differ- 
ent prescriptions have been worked out for modeling final state showers [[L2|] as opposed to 
initial state (backward) showers | I5fl . At this stage, in a full simulation, the explicit par- 
ton emissions would be combined with a hadronization model which converts the partons 
into detectable particles. Our calculation does not include hadronization. The inclusion of 
hadronization should not alter our conclusion that multiple parton emission in the initial 
state can qualitatively account for the discrepancy between theory and experiment in direct 
photon production. 

While the PS prescription for LL calculations is straightforward, the prescription for 
merging PS with NLO calculations is not. One problem is that the shower emission from a 
2 — > 2 subprocess may be double counted by the exact emission of an extra parton in the 
2 — ^ 3 subprocess. Another problem is that, to be consistent, NLO dynamics should be used 
to govern the parton shower development. We use initial and final state shower algorithms 
consistent with LL dynamics, although we use the NLO parton distribution functions in 
our calculation of initial state shower probabilities. Consequently, our calculation is not 
consistent to NLO: the PS algorithm here should be regarded only as a parametrization of a 
fully consistent NLO PS program. From a practical standpoint, the error induced by using 
only collinear dynamics in the PS algorithm in the first place should be far larger than the 
error induced by neglecting NLO corrections to the underlying collinear shower dynamics. 
Our goal here is to demonstrate that multiple parton emissions may be responsible for the 
discrepancy between data and theory at low transverse momentum. 

To avoid the double counting problem, we restrict shower development to the 2 — > 3 
subprocesses in which all momentum vectors are large and well separated. One can view 
a Monte Carlo NLO calculation as a sort of truncated parton shower, with only a single 
extra parton emission, but which is performed exactly to 0{aa 2 s ). In this case, the 2^2 
contributions, which include various soft and collinear terms for which the starting shower 
virtuality would be tiny, would never shower. If the starting shower virtuality is appropri- 
ately chosen for the 2^3 subprocesses, then only energetic, well-separated 3-body final 
states will develop a parton shower. Thus, the third parton of the 2 — > 3 subprocess can be 
viewed as the first of the potentially multiple emissions, but which is performed using exact 
instead of collinear dynamics. 

In our calculation of direct photon production, we have started with the NLO calculation 
of Ref. |7[ merged with the PS along the lines of the preceding discussion. Our computer 
program generates 2 — ^ 2 subprocesses, which frequently have negative weights, along with 
2 — > n processes, with positive definite weights, but where n > 3. Crucial to our calculation 
is the stipulation of the starting virtualities for the parton shower. 

A naive choice of starting virtuality t v , such as \t v \ = ^^(7), (with n ~ 1) does not 
ensure that the 3-parton final state is well separated. This choice leads to large amounts of 
showering even for soft or collinear configurations. One example of allowed showering with 
I it, I = np\ is a highp^ photon recoiling against two nearly collinear partons, with |t 45 | > 8 C S\2 
but still small. This is a region of phase space where the 2 — > 2 and 2^3 contributions 
at a specific Pt{i) niay cancel. Since showering is implemented only in the 2^3 processes 
and may result in a boosted Pt(i) for the 2 — > 3 contribution, the required cancellation may 
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not occur. This introduces a dependence on S c (and S s for other configurations) which is 
unphysical. In our procedure for merging NLO with PS, we minimize (but never completely 
eliminate) the dependence of results on variations of parameters. 

To minimize the dependence of results on S s and 5 C , we set the starting virtuality for 
initial state partons to \t v \ = c^min(|tjj|, s^) for i,j = 1 — 5, namely, the minimum of all 
invariants formed by the five momenta in the 2^3 process, up to a multiplicative constant 
c v . With this prescription, any nearly soft or collinear emissions in the 2 — > 3 subprocess 
will result in small starting virtualities, and a small probability to shower. Only energetic, 
well separated 2^3 subprocesses will develop a significant parton shower in the initial 
state. The final state showers are initiated with starting virtuality S12. Final state showers 
do not change Pt(j) relative to the unshowered calculation; they can, however, affect the 
number of final state photons passing the isolation cut. 



III. CALCULATIONAL RESULTS AND COMPARISON WITH DATA 

Direct photon production data from a variety of fixed target and collider experiments 
have been tabulated as a function of xt(j) = 2pT( r y)/y/s in two recent studies pl,|10 



To compare against NLO calculations, it has proven convenient to plot the quantity 
(Data — Theory) /Theory. Thus, data in perfect agreement with theory would lie along 
the y = horizontal line. In Ref. |§, a common trend amongst the various experimental 
data sets was noticed, when compared against NLO QCD. For almost all data sets tabu- 
lated, the low xt(j) range was underestimated by the theory (NLO QCD). In Ref. fL0| , 



the authors were able to improve somewhat the data vs. theory discrepancy by adjusting 
independently the factorization and renormalization scales. Nevertheless, the discrepancy 
between data and theory persists. 

In Fig. la, we show (Data — NLO)/NLO vs. Xt(t) f° r data from the CDF experiment 
16| at the Fermilab Tevatron, using pp collisions at yfs = 1.8 TeV. The data points are 



taken from Ref. ||, where the NLO distributions are calculated using the CTEQ2M PDF's 



17j evaluated at the renormalization/factorization scale /i = Pt(i)- The large enhancement 
of data over theory can be seen below Xt(j) ~ 0.05, which corresponds to Pt(i) ~ 45 GeV 
at the Tevatron. Our calculation employs the same scale choices as Ref. ||, but updated 
CTEQ3M PDF's In keeping with CDF cuts, we require the photon pseudorapidity 

(77(7) I < 0.9, and a photon isolation cut which requires that the sum of energy, projected 



transverse to the beam axis, (E^) of parton i within a cone of size AR = J (Ar/) 2 + (A0) 2 = 
0.7 satisfy 

Y,E l T < 2 GeV. 

i AR=0.7 

These two cuts are also used in Figs. 2 and 3 below. 

To minimize differences due to parton distribution choices, etc., rather than comparing 
the data to our NLO calculation merged with parton showers (NLO©PS), we show the effect 
of showering as an excess or deficit relative to the unshowered NLO calculation. In Figs, lb 
and lc, we show the relative a>r(7) distributions (NLO©PS-NLO)/NLO where the initial 
state virtuality is chosen with c v = 4. In our calculation, we have run for subprocess photon 
Pt{i) > 4 GeV, since the matrix elements are singular as Pt(i) 0; the results do not 
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change noticeably if instead we use Pt(j) > 2 GeV. Fig. 16 employs 5 S = 105 c = 0.1, and 
Fig. lc has S s = 10S C = 0.02. We see in Figs. 16 and lc that the incorporation of the PS has 
led to an enhancement of the relative xt(t) distributions at Xt(j) ~ 0.02 of about 30 — 40%, 
and hence is in accord with the data for the low range of £t(7)- The enhancement has been 
traced to the fact that a small fraction of the large population of very low £t(7) photons gets 
boosted up to higher energies by recoiling against the multiple parton emissions. Although 
the enhancement at low 0^(7) from the NLO©PS calculation is similar for the two cases, 
the large relative xt(i) distributions show a deficit of 10 — 20%. The high xt(j) deficit is 
due to the effect of the photon isolation cut. 

For very high energy events, there can still exist significant shower virtualities for events 
with quasi-soft or collinear partons, which introduces a slight dependence on 5 S and 5 C . 
There is some enhancement in showering for very high energy events, which leads to fewer 
isolated photons, and a net diminution of signal due to the isolation cut. 

If we modify the initial shower virtuality magnitude by varying c v , we find that a choice 
of c v ~ 1 results in modest enhancements of the low xt(7) region by only ~ 10%. Choosing 
c v as high as c v ~ 9 yields enhancements typically around 80%. Also, we have investigated 
how the results change by changing the initial state shower cutoff virtuality choice from 
^min = 5 GeV 2 to t min = 3 GeV 2 . The latter variation yields typically a 20% effect. In 
spite of these various uncertainties, the overall qualitative trend of enhanced cross section 
at xt{^) 0.06 persists in all the cases we have examined. 

In Ref. P, it was noted that an ad-hoc Gaussian smearing of the subprocess px led 
to improved agreement between theory and data. In Fig. Id, we additionally introduce 
Gaussian smearing (GS) to both 2 — > 2 and 2 — > 3 processes, with average transverse 
momentum zero and width a = 1 GeV. The overall enhancement of the NLO©PS at £t(7) ~ 
0.02 remains, but with some slight additional enhancement for NLO©PS©GS at even lower 
xt(j) values. The small effect of the Gaussian smearing at CDF is not surprising since the 
average boost generated by the PS algorithm is ~ 2.5 GeV. 

In Fig. 2 a, we show data from the UA2 experiment [pl| (pp collisions at a/s = 630 
GeV) compared with NLO QCD, for scale choice /1 = pt{ , j)/2. Here we use a photon pt 
cutoff of Pt{i) > 2 GeV. Again, we see that data exceeds theory by ~ 40%, although 
this time for xt(j) ~ 0.05 (corresponding to Pt(j) ~ 16 GeV). In Fig. 26, we plot the 
NLO©PS result, using again the initial state virtuality choice c v = 4, and for /i = p T ( r y)/2 
and 5 S = 105 c = 0.02. Our merged NLO©PS calculation gives an enhancement of ~ 20% 
above NLO results for xt(j) ~ 0.05. Although the CDF and UA2 calculations start with 
similar virtualities, the relatively higher value of Feynman-x in the UA2 case leads to lesser 
amounts of initial state PS radiation. This can be offset to some extent by choosing a higher 
starting virtuality, c v = 9, shown in Fig. 2c. The increase in virtuality leads to a rise in our 
calculation to about 40% above NLO expectations, in accord with the data. Finally, in Fig. 
2d, we include as well the Gaussian smearing, which leads to some additional enhancement 
at low xt(j)- 

Finally, we turn to much lower energy pp collider results from experiments at the CERN 
ISR at t/s = 63 GeV. In Fig. 3a, we show the data from the R806 experiment compared 
with NLO QCD for /i = p T (j)/2. Using the same scale /1, and including parton showers, 
we show in Figs. 36 and 3c, the comparison (NLOffiPS-NLO)/NLO for c v = 4 and c v = 9, 
respectively. We have lowered the photon px cutoff here to Pr(l) > 1 GeV. Because of the 
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large values of parton x and small virtualities, at this energy, there is very little showering, so 
that perturbative multiple parton emission as described by the PS algorithm cannot explain 
the data/theory discrepancy. However, Gaussian smearing on the order of 1 GeV can be a 
large effect at this energy, where Xt(j) = 0.1 — 0.4 corresponds to Pt(i) = 3 — 13 GeV. In 
Fig. 3d, we invoke as usual the a ~ 1 GeV Gaussian smearing of the subprocess transverse 
momentum. In this case, the smearing can move the low xt(j) theoretical prediction into 
rough agreement with the data. 

IV. SUMMARY AND CONCLUSIONS 

In summary, we have investigated the effects of multiple parton emissions on direct 
photon production in hadronic collisions by merging the PS technique with NLO QCD. 
For experiments at very high energy (e.g. UA2 and CDF), the extra k? smearing of the 
hard scattering subprocess induced by the multiple parton emissions can cause some of the 
relatively numerous low pt photons from NLO QCD to be boosted to higher pt values. 
Such an effect causes a shift in the predicted £7(7) distribution, thereby improving the 
agreement between theory and experiment. Our results cannot be interpreted as a QCD 
prediction due to the many uncertainties in the PS algorithm, and in our merging procedure. 
Amongst these uncertainties are the nature of the PS algorithm itself, and the prescription 
for initial and cutoff virtualities in the PS. On the other hand, our results can be interpreted 
as an existence proof that higher order effects (particularly from multiple parton emission) 
can account for the theory vs. data discrepancy. Other groups [|TQ| , [TT|1 have noted that 
the theory vs. data discrepancy can be resolved in NLO QCD mainly by using modified 
parton distribution functions. We comment that our result of an appropriately shifted sct(t) 
distribution will obtain for any choice of PDF's or hard scattering scale choices, as long as 
sufficient parton showering can be produced. Since hard scattering processes in nature are 
of course all-orders processes, one would expect at some level a discrepancy between data 
and fixed order QCD to occur. Our results show that this may already be the case for the 
direct photon Xt(j) distributions. 

For lower energy data sets (e.g. ^/s ^ 63 GeV), it is difficult to produce sufficient QCD 
radiation via the PS to improve the theory vs. data discrepancy. We do note, as in Ref. |§, 
that an intrinsic Gaussian kx smearing with width a ~ 1 GeV will push the theory in the 
right direction to match with data. Thus, the theory vs. data discrepancy can be resolved 
globally by invoking extra k? for the hard scattering partons: that k? would be primarily 
perturbative in nature for high energy data sets, but mainly non-perturbative for data sets 
taken at < 100 GeV. 
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FIGURES 



FIG. 1. x T {i) distribution for pp collisions at <Js = 1.8 TeV. We show a) (Data-NLO) /NLO 
for CDF data, b) (NLO©PS-NLO)/NLO for S s = 106 c = 0.1, c) (NLO©PS-NLO)/NLO for 
S s = 105 c = 0.02, and d) (NLO©PS©GS-NLO)/NLO for S s = 10S C = 0.02. For all plots, the 
hard scattering pt(t) > 4 GeV. 

FIG. 2. x t (t) distribution for pp collisions at ^/s = 0.63 TeV. We show a) (Data-NLO) /NLO 
for UA2 data, b) (NLO©PS-NLO)/NLO for c v = 4, c) (NLO©PS-NLO)/NLO for c v = 9, and 
d) (NLO©PS©GS-NLO)/NLO for c v = 9. For all plots, the hard scattering p T {l) > 2 GeV and 
S s = W5 C = 0.02. 

FIG. 3. xt(t) distribution for collisions at y/s = 63 GeV. We show a) (Data-NLO)/NLO for 
R806 data at the CERN ISR, b) (NLOffiPS-NLO)/NLO for c v = 4, c) (NLO0PS-NLO)/NLO for 
c v = 9, and d) (NLO0PS0GS-NLO)/NLO for c v = 9. For all plots, the hard scattering p T (^) > 1 
GeV and 5 S = 105 c = 0.1. 
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